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Аннотация 

Введение. Известно, что решение задачи теории упругости возможно в перемещениях, напряжениях, а также 
смешанным методом. Наиболее распространен метод конечных элементов (МКЭ) в форме метода перемещений, 
однако он имеет свои недостатки. Целью настоящей статьи является изучение особенностей численного решения 
плоской задачи теории упругости в напряжениях с использованием метода конечных разностей (МКР) для 
возможности его использования в качестве альтернативы МКЭ. Рассматриваются варианты выбора основной 
системы метода сил для получения граничных условий при расчете прямоугольной балки-стенки методом сеток 
и их влияние на конечный результат при вычислении напряжений. 

Материалы и методы. Расчет прямоугольной балки-стенки реализован в программе ВаШа-Я1епка, 
разработанной в ДГТУ. Выполняется сравнение решения на основе метода конечных разностей с решением 
методом конечных элементов в форме метода перемещений в программном комплексе «ЛИРА-САПР». 
Точность расчетов определяется путем выполнения статических проверок равновесия отсеченной части 
конструкции. Для рассматриваемой конструкции проанализированы пять вариантов выбора основной системы 
метода сил (ОСМС). 

Результаты исследования. Впервые установлено, что, в отличие от расчета стержневых систем, выбор ОСМС 
несущественно влияет на окончательные величины напряжений при решении плоской задачи теории 
упругости. 

Обсуждение и заключение. Сравнение результатов расчета при помощи метода конечных элементов в форме 
метода перемещений и метода конечных разностей показало некоторые преимущества последнего. Метод 
конечных разностей может выступить альтернативой методу конечных элементов в форме метода перемещений 
для решения плоской задачи теории упругости при условии автоматизации процесса построения эпюр 


внутренних усилий в контуре балки-стенки. 
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Введение. В качестве основного метода расчета строительных конструкций в настоящее время выступает 
метод конечных элементов (МКЭ) в форме метода перемещений [1-4]. 

Недостатками МКЭ в форме метода перемещений являются [5—7]: 

1. Пониженная, по сравнению с перемещениями, точность вычисления усилий (напряжений). 

2. Сложность учета разрывов значений усилий (напряжений) на межэлементных границах. 

3. Неиспользование граничных условий, выраженных в усилиях (напряжениях). 

4. «Запирание» при расчете с учетом элементов, обладающих большой сдвиговой жесткостью, когда энер- 
гия деформации от сдвига значительно превышает энергию деформации от изгиба. 

5. Содержащиеся в кинематических граничных условиях смещения конечного элемента, как жесткого це- 


лого, не связанные с работой внутренних сил. Они приводят к ухудшению обусловленности матрицы жесткости 
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всей конструкции при сгущении конечно-элементной сетки, когда необходимо повысить точность решения в 
опасных зонах или при итерационном решении нелинейных задач. 

6. Сложности при включении в конструкцию недеформируемых элементов и при расчетах геометрически 
изменяемых систем в виде кинематической цепи. 

Альтернативой МКЭ при решении плоской задачи теории упругости (ТУ) выступает метод конечных раз- 
ностей. При этом в качестве искомой функции выступает функция напряжений. 

Классическим объектом, для которого имеет место плоская задача ТУ, является балка-стенка — тонкая 
пластина, подверженная действию нагрузок, приложенных в ее срединной плоскости. 

В качестве граничных условий при расчете балки-стенки выступает рамная аналогия. Если рассмотреть 
контур балки-стенки как раму, то функция напряжений на контуре будет равна изгибающему моменту, а ее 


производная по нормали — продольной силе: 


ф=М; 

дф 

— =М. 1 
> (1) 


Замкнутый контур балки-стенки при рассмотрении его как рамы является трижды статически неопреде- 
лимым, однако эпюры М и М можно строить в основной системе метода сил (ОСМО). При этом в зависимости 
от выбора ОСМС эпюры будут различными, а значит и будет отличаться функция напряжений ф. Однако окон- 
чательные значения напряжений должны быть одинаковы, поскольку все функции напряжений, получаемые 
на основе различных ОСМС, отличаются на линейную функцию Ах + Ву + С. Известно, что для стержневых 
систем при решении методом сил точность результатов зависит от выбора ОСМС [8-9]. 

Целью настоящей работы выступает исследование влияния выбора ОСМС на точность определения 
напряжений в балке-стенке. В рамках поставленной цели решаются следующие задачи: 

1. Расчет балки-стенки методом конечных разностей при различных вариантах выбора ОСМС. 

2. Расчет той же конструкции при помощи метода конечных элементов в форме метода перемещений. 

3. Сопоставление результатов расчета при помощи МКР в различных ОСМС, а также МКР и МКЭ. 

Материалы и методы. Задача определения напряжений деформированного состояния балки-стенки сво- 
дится к однородному дифференциальному уравнению четвертого порядка в частных производных [10-11]: 

Г*ф = 0, (2) 
где ф — функция напряжений (функция Эри); 


9% 9* 9* 
— ———-+ —— двойной оператор Лапласа. 
дх\ дх2д22 024 д р Р 


= 
Разностная аппроксимация уравнения (2) имеет вид: 


(ба? + За + 6)ф, —- 4(1+а)(ф. +аф, +. +афа) + 
(3) 
+2а(ф. + Ф; + Фо +Ф,) + ф: + а?ф, + фи + а?фи = 0, 


где а = 42? /Дх?. 

Разностная схема (3) имеет второй порядок точности. 

Схема узлов к формуле (3) показана на рис. 1. 

Уравнение (3) составляется для всех внутренних узлов сетки. При составлении уравнений для смежных к 
контуру узлов возникает необходимость определения функции напряжений в «законтурных» точках. Для этого 
используется второе условие в (1). 

Нормальные и касательные напряжения при отсутствии объемных сил вычисляются по формулам: 

д*ф д?ф д*ф (4) 


=. == . 
922’? дх2’ дхд2 


9х 
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Или в форме метода конечных разностей: 


_ Фь- 2фк + Фа 
И 
Фа - 2Фк+Ф 
т (5) 
- _ Фо + Фе) — Фу Фь) 
т 4ДхХДЕ ° 
т 
Ф 
е а Л 
Ро 4 А Ь ® [ 
2 
Г. [6 5 
® 771 
м м] м | м] 


с 


Рис. 1. Схема узлов для составления уравнений МКР (рисунок авторов) 


В данной работе рассматривается прямоугольная балка-стенка (тонкая пластина) из изотропного материала, 
соответствующего закону Гука, загруженная нормальной к контуру поверхностной нагрузкой. При этом 
сосредоточенная сила Ё рассматривается как предельный случай равномерно распределенной нагрузки на малом 


участке балки-стенки. Толщина балки-стенки принимается равной единице. 


Требуется определить напряженное состояние пластины. Расчетная схема приведена на рис. 2. Расчет 


выполняется при а = 9,9 м; Ь = 13,5 м; 4= 30 — Р= 130 кН; Ах = 1,1 м; Д: = 0,9 м. Конструкция разбивается на 


9 интервалов по горизонтали и 15 — по вертикали. 


13 пар аа5р чер чаи чавр 119} 
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Рис. 2. Исходная схема (рисунок авторов) 
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В качестве ОСМС рассматриваются варианты, представленные на рис. 3. Опоры в узлах 113 и 152 вОСМС 
остаются нетронутыми, в вариантах (а) и (в) врезается три промежуточных шарнира в угловых точках, а в 


остальных вариантах делается разрез на одной из сторон. 


Е 
СГ. 


Рис. 3. р ОСМС и авторов) 


Расчет реализован в программе ВаКа-З{епКа, разработанной на кафедре «Сопротивление материалов» 
Донского государственного технического университета. Общее количество уравнений и неизвестных для 


рассматриваемой схемы — 112. 


11,45 
17.07 17,38 14.05 


кк << 


Р 
Рис. 4. Статическая проверка (рисунок авторов) 
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После определения напряжений выполняются две статические проверки равновесия отсеченной части. На 


рис. 4 показано равновесие верхней и правой отсеченной части при расчете на основе ОСМС варианта (а). 


Также выполняется расчет при помощи МКЭ в форме метода перемещений в программном комплексе 


«ЛИРА-САПР» при том же количестве узлов сетки. При расчете в ПК «ЛИРА» используются плоские 


прямоугольные конечные элементы балки-стенки (КЭ 21) с двумя степенями свободы в узле (линейные 


перемещения и и и вдоль осей хи 7). 


Результаты исследования. В таблицах 1-2 представлены полученные в результате расчета методом 


конечных разностей с использованием различных ОСМС значения напряжений по линиям 144—129 


(горизонтальное сечение) и 117 - 156 (вертикальное сечение). Также в указанных таблицах приведены результаты 


на основе МКЭ. 


Таблица 1 
Значения напряжений од’ и т‚„ по линии 144—129 
Вариант Номер узла 144 49 50 51 5 

0» -24,924 -37,231 -46,807 -50,975 -50,023 
а 

ха 0 -11,445 -17,068 -17,384 -14,047 
й 0, -24,946 -37,251 -46,818 -50,977 -50,019 

тг 0 -11,451 -17,076 -17,393 -14,054 

Г -24,948 -37,251 -46,817 -50,977 -50,019 
В 

Тя 0 -11,451 -17,077 -17,393 -14,054 

0, -25,17 -37,324 -46,783 -50,883 -49,918 
г 

Тоха 0 -11,501 -17,133 -17,424 -14,05 

[ей -25,051 -37,284 -46,793 -50,931 -49,977 
д 

тг 0 -11,487 -17,104 -17,399 -14,043 

[ор -30,7 -36,1 -45,4 -50,5 -50,2 

ЛИРА-САПР 
хх -6,33 -10,8 -17,1 -17,9 -14,4 
Вариант Номер узла 53 54 55 56 129 

0», -44,513 -34,072 -17,462 6,209 34,671 
а 

Тоха -8,575 -2,431 2,363 3,261 0 

0, -44,506 -34,064 -17,453 6,219 34,685 
6 

ха -8,58 -2,435 2,361 3,26 0 

0, -44,506 -34,064 -17,453 6,219 34,685 
В 

ха -8,58 -2,435 2,361 3,26 0 

Гей -44,433 -34,034 -17,471 6,152 34,559 
г 

хе -8,549 -2,392 2,402 3,287 0 

Гей -44,48 -34,056 -17,463 6,191 34,638 
д 

хе -8,561 -2,414 2,379 3,271 0 

0» -44.4 -32,9 -15 8,83 22 

ЛИРА-САПР 
хх -8,33 -1,81 2,49 2,62 1,54 
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Таблица 2 
Значения напряжений д. ит,» по линии 117-156 
Вариант Номер узла 117 4 12 20 
Ох -125,231 -99,051 -64,507 -34,823 
ы Е 0 -1.958 -5,954 -9.754 
б бх -125,212 -99,038 -64,499 -34,821 
хи 0 -1,964 -5,962 -9,763 
©х -125,213 -99,038 -64,499 -34,821 
- хи 0 -1,963 -5,962 -9,763 
| ©х -125,251 -99,067 -64,526 -34,849 
хи 0 -1,965 -5,965 -9,766 
©х -125,225 -99,049 -64,51 -34,834 
ь хи 0 -1,964 -5,963 -9,764 
бх -114 -99,8 -88,3 -37,3 
ны хи -0,582 -1,86 -5,22 -9,33 
Вариант Номер узла 28 36 44 52 
бх -14,398 -2,907 1,454 0,406 
ы хи -12,224 -13,406 -13,821 -14,047 
©х -14,398 -2,911 1,448 0,397 
о хи -12,233 -13,414 -13,829 -14,054 
©х -14,398 -2,911 1,448 0,397 
ь хи -12,232 -13,414 -13,829 -14,054 
бх -14,429 -2,942 1,422 0,386 
ь хи -12,235 -13,415 -13,827 -14,05 
0х -14,414 -2,929 1,429 0,383 
к хи -12,233 -13,413 -13,824 -14,043 
бх -14,6 -1,88 1,76 0,61 
а То -12,6 -13,9 -14,3 -14,4 
Вариант Номер узла 60 68 76 84 
ь бх -4,885 -13,717 -25,506 -39,38 
хи -14,564 -15,675 -17,342 -18,897 
б ©х -4,897 -13,732 -25,523 -39,396 
хи -14,57 -15,679 -17,342 -18,89 
бх -4,897 -13,732 -25,523 -39,396 
В та -14,57 -15,679 -17,342 [18,89 
бх -4,877 -13,667 -25,408 -39,247 
ь хи -14,568 -15,688 -17,379 -18,962 
©х -4,898 -13,706 -25,455 -39,283 
ь хи -14,551 -15,653 -17,318 -18,892 
бх -4,14 -13,4 -25,8 -40,6 
о Тк» -15 -16,4 -18,5 -20 
Вариант Номер узла 92 100 108 156 
бх -53,872 -66,502 -73,431 -70,302 
ь ты -18,924 -15,761 8.761 0 
б бх -53,881 -66,496 -73,404 -70,264 
хи -18,91 -15,741 -8,744 0 
в бх -53,881 -66,496 -73,404 -70,264 
хи -18,91 -15,741 -8,744 0 
ы бх -53,735 -66,382 -73,386 -71,111 
хи -18,983 -15,723 -8,633 0 
бх -53,759 -66,469 -73,571 -70,422 
- хи -18,968 -15,842 -8,732 0 
бх -55,4 -86,9 -71,8 -72,2 
Вы ее -18,8 -14,1 -7,1 -3,15 
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В таблицах 3—4 представлены полученные погрешности статических проверок равновесия левой и верхней 
отсеченной части при различных вариантах ОСМС, а также при расчете с помощью МКЭ в программном 


комплексе «ЛИРА-САПР». 


Таблица 3 
Погрешности статической проверки равновесия левой отсеченной части 
Погрешность ХХ РУ ХМ 
а) 0,313 % 0,478 % 0,0805 % 
6) 0,001 % 0,476 % 0,133 % 
в) 0,001 % 0,476 % 0,133 % 
г) 0,0036 % 0,461 % 0,0827 % 
д) 0,00103 % 0,465 % 0,16 % 
ЛИРА-САПР 0,0405 % 0,611 % 0,522 % 
Таблица 4 
Погрешности статической проверки равновесия верхней отсеченной части 
Погрешность ХХ РЯ ХМ 
а) 0,979 % 0,0017 % 0,454 % 
6) 0,99 % 0,73 % 0,464 % 
в) 0,99 % 0,73 % 0,464 % 
г) 1,007 % 0,0056 % 0,453 % 
д) 0,988 % 0,0017 % 0,462, % 
ЛИРА-САПР 0,942 % 0,0069 % 0,542 % 


Из таблиц 1-2 видно, что напряжения при расчете методом конечных разностей при выборе основной 
системы метода сил по схемам «а—д» отличаются незначительно. Среднее расхождение результатов, полученных 
при помощи МКР, с решением на основе метода конечных элементов составляет 13,9 % по напряжениям д„, 
7,2 % по напряжениям д, и 7,7 % по напряжениям т,„,„ (без учета крайних узлов). Следует также обратить 
внимание на то, что в лежащих на контуре балки-стенки узлах 144, 129, 117, 156 при расчете методом конечных 
элементов касательные напряжения т,„ оказываются отличными от нуля, что противоречит условиям на 


поверхности: 
Ху, = а, [+т,. ‘п, 
(6) 
И, = тих '[+0,.п, 
где Х, и 7, составляющие поверхностной нагрузки по осям х и 2; |= с05(х,У); п = с0$(2%) — направляющие 
косинусы нормали к поверхности балки-стенки "У. 

Невыполнение условий на поверхности по касательным напряжениям при расчете в ПК «ЛИРА» можно 
объяснить тем, что напряжения в данном программном комплексе вычисляются только в центрах конечных 
элементов, и выведенные значения для узлов 144, 129, 117, 156 на самом деле соответствуют точкам, смещенным 
на Дх/2 или Д2/2. 

В то же время эта особенность никак не отразилась на погрешностях статических проверок, они 
сопоставимы с результатами, полученными при помощи метода конечных разностей региона [9]. 

Обсуждения и заключения. Полученные результаты позволяют сделать следующие выводы: 

1. Выбор основной системы метода сил практически не влияет на результат по напряжениям. 

2. Точность решения при помощи метода конечных элементов сопоставима с точностью метода конечных 


разностей, но при решении методом конечных элементов при одинаковой густоте сетки число неизвестных в узле 


в 2 раза больше (перемещения 1 и №). 
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3. Недостатком МКЭ по сравнению с МКР является невозможность точного определения напряжений 
непосредственно на поверхности конструкции. 

4. МКР может служить альтернативой МКЭ в форме метода перемещений при решении плоской задачи 
теории упругости в случае автоматизации процесса построения эпюр изгибающего момента и продольной силы 


в контуре балки-стенки, как в раме. 
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